Non-gaussian statistics and the relativistic nuclear equation of state 
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We investigate possible effects of quantum power-law statistical mechanics on the relativistic nu- 
clear equation of state in the context of the Walecka quantum hadrodynamics theory. By considering 
the Kaniadakis non-Gaussian statistics, characterized by the index « (Boltzmann-Gibbs entropy is 
recovered in the limit k — > 0), we show that the scalar and vector meson fields become more intense 
due to the non-Gaussian statistical effects (k/0). From an analytical treatment, an upper bound 
on k (k < 1/4) is found. We also show that as the parameter k increases the nucleon effective 
mass diminishes and the equation of state becomes stiffer. A possible connection between phase 
transitions in nuclear matter and the K-parameter is largely discussed. 

PACS numbers: 21.65.+f; 26.60.+C; 25.75.-q 



I. INTRODUCTION 



In the theoretical treatment of the properties of nuclear 
matter, the relativistic phcnomcnological approach devel- 
oped by Walecka [l[ , the so-called quantum hadrodynam- 
ics (QHD-I), represents one of the important approaches 
to the highly nonlinear behavior of strong interactions at 
the hadronic energy scales. This model provides a ther- 
modynamically consistent theoretical framework for the 
description of bulk static properties of strong interact- 
ing many-body nuclear systems Formally, QHD-I is 
a strong-coupling rcnormalizable field theory of nucleons 
interacting via the exchange of (isoscalar) scalar (a) and 
vector(w) mesons [l|. The model has been largely used 
in calculations of nuclear matter and finite nuclei (see, 
e.g., 0] and Refs. therein). 

An important aspect worth emphasizing concerning 
the Walecka treatment is that it is based on standard 
quantum statistical relations, i.e., Fermi-Dirac distribu- 
tions. On the other hand, as is well known, some re- 
strictions to the applicability of the standard statistical 
mechanics have motivated investigations of power-law or 
non-gaussian statistics, both from theoretical and experi- 
mental viewpoints. In this concern, the Tsallis nonexten- 
sive statistical mechanics [3| and the extensive general- 
ized power-law statistics developed by Kaniadakis H are 
the most investigated frameworks. Several consequences 
(in different branches) of the former framework have been 
investigated in the literature 6], which includes systems 
of interest in high energy physics, namely, the problem 
of solar neutrino Q, the charm quark dynamics in a 
thermal quark-gluon plasma for the case of collisional 
equilibration [8j|, interpretations for central Au-Au colli- 
sions at RHIC energies in a Relativistic Diffusion Model 
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(RDM) [H, among others (see, e.g., The Kani- 

adakis non-gaussian statistics in turn is characterized by 
the K-entropy (see Sec. Ill) that emerges naturally in 
the framework of the so-called kinetic interaction prin- 
ciple Q . Several physical features of the K-distribution 
have also been theoretically investigated as, for instance, 
the self-consistent relativistic statistical theory [5], non- 
linear kinetics [l(|, and the H-theorem from a general- 
ization of the chaos molecular hypothesis [Til ]. 

In this paper, by following our previous results [l2| , we 
adopt the Kaniadakis non-gaussian statistics and study 
the effects of the K-generalization for the Fermi-Dirac 
(FD) distributions on the QHD-I theory, especially what 
concerns its effects on the equation of state (EoS) of the 
nuclear matter. We also investigate the influence of these 
non-gaussian effects on phase transitions in the nuclear 
matter, as discussed in the standard context by Theis et 
al. 0. 

This paper is organized as follows. In Sec. II, we 
present the basic formalism of the mean field theory of 
QHD-I, which is important for the calculation of some 
quantities of nuclear matter. A brief review of Kani- 
adakis statistics is presented in Sec. III. In Sec. IV the 
convergence of the calculation is discussed and it is shown 
that the allowed values of k lie in the range < n < 0.25. 
Our main results are discussed in Sec. V. In Sec. VI, we 
consider the connection between the k— generalisation of 
the QHD-I theory and the phase transition in nuclear 
matter, as discussed in the paper (l3| . We summarize 
our main conclusions in Sec. VII. 



II. BASICS OF QHD-I 

The Lagrangian density describing the nuclear matter 
reads [l[ 

C = i>[(n^ - g w ^) - (M - g a a)]^ 
+ \{d^a - mW) - \u^ v + \ml^ , (1) 
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which represents nuclear matter composed by nucleons 
coupled to two mesons, namely, the a and u> mesons (for 
details see reference [l|). 

Applying standard techniques from field theory and 
the mean-field approach, we obtain the scalar density 
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where M* is the effective mass 



[n(f,T) + n(v,T)]d 3 k , (2) 



M* =M 



j0 - = M - ^p s 



(3) 



The baryon number density, the energy density and pres- 
sure are given, respectively, by 



Qb = 



7n 
(2tt) 3 



[n(v,T) - n(v, T)}d 3 k, 
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and 
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i{v,T) = 



[n(is,T) + n(v,T)}d 3 k , (6) 
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,J3[E'(k)- v \ , j 



(7) 



^ T ) = ^MTT- (8) 

are the usual FD distributions for baryons and anti- 
baryons, with E*(k) = \/k 2 +M* 2 , (3 = l/k B T. The 
parameter v = p — g u u)Q = p — [g^ j m u} ')' 1 Qb is the effec- 
tive chemical potential, and 7^ is the multiplicity factor 
(77V = 2 for pure neutron matter and 7^ = 4 for nuclear 
matter) . 

Additionally, to obtain the results described in Sec. V, 
we use for the coupling constants the values of reference 
[l|, namely 1 , 

= 1 1 .798 fm 2 and ( ^] = 8.653 fm 2 , (9) 



which are fixed to give the bind energy -Ebind = —15.75 
MeV and k F = 1.42 fm -1 . In Sec. VI other values of 
(<7o-/to ct ) 2 are considered. 



1 For the purpose of the present work, the values given in Eq. 
||9} suffices to investigate the effects of the power-law statistics 
in neutron and nuclear matter. Variations of the coupling con- 
stants, within the acceptable values given in current literature, 
do not qualitatively affect the conclusions. 



III. NON-GAUSSIAN FRAMEWORK 

Recent studies on the kinetic foundations of the so- 
called K-statistics led to a power-law distribution func- 
tion and a K-entropy which emerges naturally in the 
framework of the kinetic interaction principle (see, e.g., 
Ref. @). Formally, the K-framcwork is based on the 
K-exponential and the K-logarithm functions which are 
defined as M 



exp K (/) = ( V / l + « 2 / 2 + «/) 1/K , 



m/) = (r - r K )/2*, 



(10) 



(11) 



ln K (exp K (/)) = expjln re (/)) = /. (12) 

The K-entropy associated with this K-framework is given 
by 



S K (f) = - J d 3 pfln K f, 



(13) 



which fully recovers standard Boltzmann-Gibbs entropy, 
S K =o(f) = ~ J f In f d3 P, in the limit «->0. 



A. Quantum Statistics 

We recall the main aspects of the connections be- 
tween the quantum statistics and Kaniadakis framework. 
Specifically, the main result is that, for values of the k 
index lying in the interval [—1,1], a K-generalized quan- 
tum distributions for fermions and bosons can be written 
as y| 



n K (fi,T) 



e K (/3(e-t*))±l' 



where e K reads 

e K (x) = ( \/l + n 2 x 2 + Kx^j 



l/re 



(14) 



(15) 



and x = /3(e — p). In particular, we observe that, 
e_ K = e K , and that in the K — * limit, the standard 
FD distribution, n(p,T), is recovered. As physically ex- 
pected, as T —>■ 0, n K (p,T) — > n{p,T). This amounts 
to saying that for studies of the interior of neutron stars 
(where, in nuclear scale, T ~ 0) we do not expect any 
power-law statistic signature. On the other hand, in 
heavy ions collision experiments or in the interior of pro- 
toncutron stars, with typical stellar temperatures of sev- 
eral tens of MeV (1 MeV= 1.1065 x 10 10 K), power-law 
statistics effects can be relevant. In order to study such 
effects, in the next Section we combine Eqs.Q-© with 
the generalized FD distributions given by Eqs. (fl~4"]) and 

G5). 
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FIG. 1: The effective nucieon mass and the vector meson fieid of pure neutron matter (7jv = 2) as function of temperature for 
different values of the parameter k. Panel (a): the self-consistent nucieon mass at vanishing baryon density. Panel (b): the 
vector meson field at nonzero baryon density corresponding to v — 540 MeV. 



IV. ^-STATISTICS AND QHD-I 

Since the above function e K (x) is a deformed exponen- 
tial function, we must verify the mathematical conver- 
gence of the integrals in Eqs. ([2])-© when considering 
the K-distribution (| 14|) . Now, taking into account that 



(16) 



e K (x) — > [2kx) 1 ' k , 



in the limit x » 1, we obtain the asymptotic behaviour 
for the integrals appearing in Eqs. AS])-®, i.e., 



M* 



M* d 3 k 



E*(k){e K [E*(k) ± v] + 1} (2k/3)V« fc 1 /* ' 

(17) 



PB 



d 3 k 



e K [E*{k) -u] + \ (2k/3)V k feV« 



E*{k) d 3 k 



e K [E*(k)-v]+l 



k 2 d 3 k 



fc 4 



(18) 



(19) 



E*(k){e K [E*(k) 



1]} 



(2k/3)V« feV«' 
(20) 

From Eqs. (pT7)> - (|20[) the general asymptotic behaviour 
can be summarized by k N /k}l K , for N — 2,3, 4. In order 
to have k N /fc 1 / K — > 0, when k — > oo, we find that 
1/k > N, from which we obtain 



At high temperatures (T — > oo), the analytic solution 
to Eq. ([3]) can be written as 



M* 



M 



(22) 



where is given by Eq. (|24|) . 

Several limiting cases of the EoS are of interest: 

1. The baryon distribution becomes a step function 
n K (u,0) = 8(k F - |k|) in the limit T —> for any 
value of pb- 

2. The system becomes degenerate in the limit ps — > 
oo at any T. 

3. As T — > oo for any vaue of /Ob, an equation of state 
similar to that of a black body is obtained: 



7jv 



where, in Eqs. (22) and (|2"3"j) . 

Jo e K (z) + 1 



(23) 



(24) 



When K^Owe have that £i(/t — > 0) — > 7r 2 /12 and 
£j(re — > 0) — > 77r 4 /120, recovering the limits of Ref. 



V. THE PHASE STRUCTURE AT NONZERO 
BARYON DENSITY 



K < 



N 



(21) 



Note that, to simultaneously satisfy the convergence of 
all integrals in Eqs.©-©, we find that < n < 1/4. 



EoS 



In order to study the effects of the non-gaussian frame- 
work on QHD-I theory, we have calculated numerically, 
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FIG. 2: Isotherms of nuclear matter (7^ = 4) equation of state at finite temperatures for different values of the parameter k. 
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FIG. 3: Panel (a): Isotherms of nuclear matter (7^ = 4) at T = 15 MeV for several values of the parameter k. Panel (b): 
Isotherms of nuclear matter corresponding to the critical temperatures and parameters k c given in the data. The curves and 
the data are in the same order from top to bottom. Panel (c): the critical parameter n c as funtion of the critical temperature 
in the range 18.8 MeV < T c < 20.2 MeV (including the data of Panel (b)). 



for several values of temperatures and of the parameter k, 
the effective nucleon mass, the vector and scalar mesons 
fields for pure neutron matter (jn = 2), as well as the 
EoS for nuclear matter (7JV = 4). 

Fig. [1] shows the non-gaussian effects on the effec- 
tive mass and vector meson for a pure neutron mat- 
ter. In Panel (a), for qb = v = 0, we note that the 
higher the parameter k, the smaller the effective mass 
M* (and, consequently, the higher the scalar field a, since 
g a <7 = M — M*). Such an effect may be physically un- 
derstood in that at a given temperature, the scalar den- 
sity, as source for the scalar mesons, increases with the 
increasing of the parameter n. Thus, the attraction of 
the nucleons, mediated by the scalar mesons, becomes 
stronger, reducing the effective mass. The behaviour of 
vector meson field g^a — {g^/ra^) 2 qb is shown in Panel 
(b). 

The non-gaussian effects on the EoS of symmetric nu- 
clear matter are shown in Fig. [21 The panels display the 
logp — loge plane for selected values of k. The results 
are ploted for arbitrarily chosen values of temperature, 



T =18 MeV, 20 MeV, and 30 MeV. In reality, the motiva- 
tion for this choice is of astrophysical interest, e.g., in the 
study of protoneutron stars. Clearly, the non-gaussian ef- 
fect is manifested in the increasing of the pressure with 
the values of k making the EoS stiffer. 

B. Phase transitions 

Another effect of the power-law statistics on nuclear 
EoS of QHD-I concerns the phase transitions. From Pan- 
els (a) and (b) of Fig. 3, we see that the first order phase 
transition may be eliminated by the variation of the pa- 
rameter k. This fact can be easily visualized in Fig. [3ja), 
where isotherms at T = 15 MeV are plotted for several 
values of k S [0, 1/4). Note that, for increasing values 
of K, the dip in the region of thermodynamical instabil- 
ity becomes smaller, vanishing at the turning point that 
defines the critical values of thermodynamical quantities 
(T c , p c , etc.). 

We also note that for our choice T = 15 MeV the 
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FIG. 4: The solution M* /M of equation of ([3]) for nuclear matter (7 = 4) at vanishing baryon density as a function of 
temperature for the same value of C| and two different values of the parameter k. The right Panel shows the same results but 
with a stretched temperature region around the point of phase transition. 



upper value of n (near the 1/4 limit discussed earlier) 
is not sufficient to eliminate the first order phase tran- 
sition [16]. On the other hand, first order phase tran- 
sitions can be eliminated for values of T in the region 
18 MeV < T < 20.2 MeV. In this interval, all tem- 
peratures can be made critical. This amounts to say- 
ing that a (critical) parameter k c can be determined 
in order to yield a turning point in the isotherm at a 
given temperature. This is illustrated in Fig. [3]Jb) where 
several isotherms are displayed. InFig. [3jc) values of 
plotted as function of temperature in the range 
(18.8 MeV < T < 20.2 MeV). 



VI. THE PHASE STRUCTURE AT ZERO 
BARYON DENSITY 

Ref. [TH discusses experimentally a reproduction of the 
observed binding energy and density of nuclear matter in 
an area characterized by a line in the coupling-constant 
plane, where part of this line defines a system exhibiting 
a phase transition around T c ~ 200 MeV. As matter of 
fact, a different sets of coupling constants in the coupling- 
constant plane were considered 2 . 

In this Section, by considering the same arguments of 
Ref. [13l ] , we explore the non-gaussian phase structure of 
the effective Lagrangian at vanishing chemical potential 
and baryon density (qb = v = Q). To this end, we first 
consider the range of values for the coupling constant C| 
given by 

C% = {M^-fM 2 (25) 



in the coupling-constant plane shown in Fig.(l) of 
Ref. [l3| • At vanishing chemical potential the terms with 
the baryon density do not appear in Eqs. © and ©. In 
what follows, we have taken as an example Cl = 365. 

Theis et al. [HI showed that the order of phase tran- 
sition is strongly dependent on the actual value of the 
coupling constant Cf . However, in Sec|V] it is shown 
that we can avoid Maxwell construction by the variation 
of the parameter k. Since the order of transition depends 
on Cg and K, the natural question is whether there exists 
some relation between C% and k at zero baryon density. 

In Fig. [H the sudden drop in M* around T ~ 185 
MeV determines the abrupt rise of the energy density 
and pressure. We note that for k = (Fermi-Dirac 
statistics) the self-consistency equation ^ has three so- 
lutions around T ~ 185 MeV imposing a sudden rise in 
the energy density and a peak in the specific heat. This 
behavior is shown in Figs. ([5])- 0, where the tempera- 
ture dependence with total energy density, pressure and 
specific heat divided by the corresponding high Stephan- 
Boltzmann temperature limit (for k = 0), given by 
Eq. (f2"3")l . are displayed. The behavior of the curves for 
k = characterizes a first order phase transition with the 
pressure curve crossing itself twice at T ~ 185 MeV. We 
observe that in this region the energy density and presure 
are also triple valued and that the specific heat is neg- 
ative. The value of k = 0.0855 was obtained by requir- 
ing that energy density and pressure to be single valued 
characterizing a second order phase transition with non- 
negative specific heat. This allows to obtain a relation 
between Cj and explained below. 

A. C§ — « relation 



2 For a different sets of coupling constant, the mean field solutions 
provide the nuclear binding energy —16 < e < —15 MeV at 
equilibrium densities 0.14 < p eq < 0.19 fm — 3 . 



In order to obtain such a relation, we investigate the 
thermodynamical behavior of the nuclear matter for sev- 
eral values of the coupling constant Cl in the coupling- 
constant plane as follows. For each value of C|, the cor- 
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FIG. 5: The total energy density et and the scalar-field energy density e a divided by the k — Stefan-Boltzmann limit ej 
as a function of temperature, at zero baryon density of nuclear matter (7 = 4). The same value of C| is considered for two 
different values of the parameter k. In the right Panel the same results in the stretched temperature region near the transition 
point. 
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FIG. 6: For nuclear matter (7 = 4) at vanishing baryon density, the total pressure divided by the corresponding k = Stefan- 
Bolztmann limit Pf B as function of temperature. The same value of Cf is considered for two different values of the parameter 
k. The right panel shows the first order phase transition point for k — (solid curve). For k = 0.0855 (dotted curve) the phase 
transition is of second order. 



responding parameter k, for which the transition is of 
second order, is determined. Thus, a curve in the n x Cf 
plane is obtained as shown in Fig. [5[a) ■ Below this curve 
the phase transitions are of the first order and above it 
the thermodynamical behavior is smooth. Let us now 
show how the calculation is made via specific heat. 

Differently from the treatment discussed in Ref. [l3|, 
the mathematical structure of the self-consistency equa- 
tion in our approach is not simple, so that the calculation 
must be done numerically. We observe that the specific 
heat calculated from Eq. ([5]) is linear in dM* /dT. So, 
whenever there is a sudden fall in M*, we see a peak in 
the specific heat. By writing 



we can see from Eq. ([3]) that 

dT ( „oo n l k K} k 2 dk poo n*" 1 M* 2 dk 

1 + ZC M ' I Jo B »(fe) Jo E»(k) 

(27) 

where Cm* — (9a/m a ) 2 ^ N /-K 2 = Cg/M 2 . The singu- 
larities of dM* /dT lie in the curve determined by the 
vanishing of the denominator. Using Eq. ([3]), this condi- 
tion becomes 



M - 2C M , I nl K} 



dk 



Vk 2 + M* 2 



, (28) 



where 



H 



de 
dT 



de dM* 
dM* dT 



(26) 



,{«} 



e K {/3(Vk 2 + M* 2 -v)] + l 



(29) 
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FIG. 7: The specific heat of nuciear matter (7 = 4) at zero baryon density divided by the corresponding k = Stefan-Boitzmann 
limit as function of temperature for the same value of C% and two different values of the parameter k. In the right Panel, the 
same results in the stretched region around the phase transition point. For k = and n = 0.0855 the phase transitions are, 
respectively, of the first and second order. 



Note that, for n = 0, we fully recover Eq. (18) of Ref. p^ |. 
The number of intersections of the solutions obtained 
from Eq. I|28p and the self-consistency equation given by 
Eq. determines how the decoupling happens. We have 
a first or a second order phase transition, respectively, for 
two or one intersections. If there is no intersections, the 
thermodynamical behavior is continuous. The numerical 
results shown in Fig. [5Ja) can be summarized as follows: 

1 . For k lying below the k x C| curve the phase tran- 
sitions are of the first order. 

2. For k lying on the k x C| curve the phase transitions 
are of the second order. The corresponding values 
of temperatures and effective masses are shown in 
Figs. Mb) and He). 

3. For k lying above the k x C s curve the decoupling 
is continuous. 

Thus, the order of phase transitions depends not only on 
the actual values of C\ but also on the values of k. 



VII. FINAL REMARKS 

In this paper, we have investigated the effects of the 
non-gaussian Kaniadakis framework on the mean field 
theory of Walecka (QHD-I) We have used, instead 
of the standard Fermi-Dirac nucleon and antinuclcon dis- 
tribution functions, the K-quantum distribution obtained 
by Kaniadakis in the framework of the kinetic interaction 
principle in Ref. []|. 

We have considered pure neutron and nuclear mat- 
ter at nonzero and zero baryon densities. In the first 



case, the non-gaussian effects on nuclear and pure neu- 
tron matter, for a considerable range of temperature, is 
to make the equation of state stiffer and to increase the 
intensity of the vector and scalar meson fields, with a 
consequent lowering of the nucleon effective mass (for in- 
creasing values of the parameter k). We believe that it 
may have consequencies in astrophysical studies, mainly 
in what concerns the calculation of masses of compact 
objects, such as protoneutron stars. 

Another interesting feature of the k-QHD-I is that, at 
temperatures in the range 18 MeV < T < 20 MeV, phase 
transitions of first order can be avoided by a convenient 
variation of the parameter k, which allows the determi- 
nation of a critical k c parameter at the turning point of 
an isotherm at a given T. 

In the second case, we have examined the phase struc- 
ture of nuclear matter at high temperature and at zero 
baryon density. The effective Lagrangian of QHD-I the- 
ory is considered for the same set of values of the coupling 
constants in the coupling constant-plane of Ref. [l3|. 
Given that, at vanishing baryon density, the only nonzero 
coupling constant is that of the scalar field, a relation be- 
tween the coupling constant C s and the parameter k is 
obtained. This relation determines, in the k x C s plane, 
regions of different thermodynamical behaviors. 
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